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Abstract. We review the spectral properties of stochastic backgrounds of 
astrophysical origin and discuss how they may differ from the primordial contribution 
by their statistical properties. We show that stochastic searches with the next 
generation of terrestrial interferometers could put interesting constrains on the physical 
properties of astrophysical populations, such as the ellipticity and magnetic field of 
magnetars, or the coalescence rate of compact binaries. 



1. Introduction 

According to various cosmological scenarios, we are bathed in a stochastic background 
of gravitational waves, memory of the first instant of the Universe, up to the limits of 
the Plank era and the Big Bang [1, 2]. In addition to the cosmological background 
(CGB), an astrophysical contribution (AGB) may have resulted from the superposition 
of a large number of unresolved sources since the beginning of stellar activity, which can 
be either short lived burst sources, such as core collapses to neutron stars [3, 4, 5, 6] 
or black holes [7, 8, 9, 10, 11], phase transition or oscillation modes in young neutron 
stars [12, 13, 14], the final stage of compact binary mergers [15, 16], or periodic long 
lived sources, typically pulsars [17, 18], the early inspiral phase of compact binaries 
[19, 20, 21, 22] or captures by supermassive black holes [23, 24], whose frequency is 
expected to evolve very slowly compared to the observation time. 

The astrophysical contribution is important for at least two reasons: on one hand, 
it carries information about the star formation history, the mass range of neutron star 
or black hole progenitors, the statistical properties of populations of compact objects 
like the ellipticity or the magnetic field of neutron stars, or the rate of compact binary 
mergers. On the other hand, it can be a foreground for the cosmological background, 
and it has to be modeled accurately to define the best frequency windows where to 
search for the cosmological background. The rest of this article will be organized as 
follows: in Sec. 2, we review the spectral properties of the AGBs; in Sec. 3, we discuss 
how their statistical properties may differ from those of the cosmological background; 
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in Sec. 4 we present two simple models, the stochastic background from magnetars 
and double neutron star coalescences, susceptible of producing a continuous stochastic 
background in the frequency band of terrestrial interferometers; in Sec. 5 we discuss 
the constraints one can expect to put on these models and on the source populations 
with the next generations of detectors. Finally, in Sec. 6 we summarize the results and 
present our conclusions. 



2. The spectral properties 

The spectrum of the gravitational stochastic background is usually characterized by the 
dimensionless parameter: 

n 9 M = -p^ (i) 

where p gw is the gravitational energy density, v Q the frequency in the observer frame and 

3H 2 

p c = the critical energy density needed to close the Universe today. For a stochastic 
background of astrophysical origin, the energy density is given by: 

tigwifo) = -\voF Vo (Vo) (2) 

where the integrated flux at the observed frequency u a is defined as: 

F Vo tyo) = j fv ("o,z) d ^dz (3) 
The spectral properties of a single source located at z are given by the fluence: 



1 dE nm , x 1 dE n 



where di = r(l + z) is the distance luminosity, r the proper distance, which depends 
on the adopted cosmology, the gravitational spectral energy emitted by a single 
source and v = v (\ + z) the frequency in the source frame. 

When the gravitational emission occurs shortly after the birth of the progenitor, 
the event rate can be associated with the progenitor formation rate: 

dR(z) _ R*(z)dV 

~dT ~ X YV~zltz~ {z) (5) 

where R*(z) is the cosmic star formation rate per comoving volume (SFR) expressed 
in M Mpc~ 3 yr _1 and A the mass fraction converted into the progenitors, assumed 
to be the same at all redshifts. The (1 + z) term in the denominator corrects the 
cosmic star formation rate by the time dilatation due to the cosmic expansion. In 
our calculations, we consider the 737 cosmology [25], with Vl m = 0.3, f2 A = 0.7 and 
Hubble parameter H Q = 70 km s _1 Mpc -1 (or equivalently h = 0.7), corresponding to 
the so-called concordant model derived from observations of distant type la supernovae 
[26] and the power spectra of the cosmic microwave background fluctuations [27]. We 
use the recent SFR of [28] derived from new measurements of the galaxy luminosity 
function in the UV (SDSS, GALEX, COMB017) and FIR wavelengths (Spitzer Space 
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Telescope), which allowed to refine the previous models up to redshift z ~ 6, with very 
tight constraints at redshifts z < 1 (and quite accurate up to z ~ 2), and normalized by 
the Super Kamiokande limit on the electron antineutrino flux from past core-collapse 
supernovas. A parametric fit is given by 



„ , , , 0.017 + 0.13,2 

R * {Z) = N + (,/3.3P (6) 
for an initial mass function (IMF) of the form (modified Salpeter A): 

°" I (li!)" 2 - 35 for rn <m< 100 {} 
with a turnover below m = 0.5 M and normalized within the mass interval 0.1 — 100 
M Q such as / mC,(m)dm = 1. The element of comoving volume in eq. 5 is given by: 

dV , . , c r(z) 2 
where E(Q,z) = ^JQa + ^m(l + z) z . Combining the expressions above, one obtains: 

n ^) = ^M (i + ^n.z) du dz (9) 

or replacing the constants G (gravitational constant) and c (speed of light) by their 
usual values: 

IW,J = " x 10-«(^) 2 A,„ /-"' «.M (10) 

/i "'o (1 + 2;) 2 E(fi, z) di/ 

where i?* is given for ho = 0.7. The upper limit of the integral, which depends on both 
the maximal emission frequency in the source frame and the maximal redshift of the 
model of star formation history (usually z nmx = 5 — 6), is given by: 

{y \f 1/ <~ ^ max 

- 1 otherwise 

Consequently, the shape of the spectrum of any astrophysical background is 
characterized by a cutoff at the maximal emission frequency and a maximum at a 
frequency which depends on the shape of both the SFR and the spectral energy density. 



3. The detection regimes 

Besides the spectral properties, it is important to study the nature of the background. 
In the case of burst sources the integrated signal received at z = from sources up to 
a redshift z, would show very different statistical behavior depending on whether the 
duty cycle: 

D(z) = J o Z f(l + z>)^(z>)dz> (12) 

defined as the ratio, in the observer frame, of the typical duration of a single event f, 
to the average time interval between successive events, is smaller or larger than unity. 
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When the number of sources is large enough for the time interval between events to be 
small compared to the duration of a single event (D » 1), the waveforms overlap to 
produce a continuous background. Due to the central limit theorem, such backgrounds 
obey the Gaussian statistic and are completely determined by their spectral properties. 
They could be detected by data analysis methods in the frequency domain such as the 
cross correlation statistic presented in the next section [29]. On the other hand, when 
the number of sources is small enough for the time interval between events to be long 
compared to the duration of a single event (D « 1), the sources are resolved and 
may be detected by data analysis techniques in the time domain (or the time frequency 
domain) such as match filtering [30, 31]. An interesting intermediate case arises when 
the time interval between events is of the same order of the duration of a single event. 
These signals, which sound like crackling popcorn, are known as "popcorn noise". The 
waveforms may overlap but the statistic is not Gaussian anymore so that the amplitude 
on the detector at a given time is unpredictable. For such signals, data analysis strategies 
remain to be investigated [32, 33], since the time dependence is important and data 
analysis techniques in the frequency domain, such as the cross correlation statistic, are 
not adapted. 



4. Models of astrophysical stochastic backgrounds 

In this section, we investigate two processes able to produce a continuous stochastic 
backgrounds in the frequency band of terrestrial interferometers. 



4-1. magnetars 

Rotating neutron stars (NSs) with a triaxial shape may have a time varying quadrupole 
moment and hence radiate GWs at twice the rotational frequency. The total spectral 
gravitational energy emitted by a NS born with a rotational period Po, and which 
decelerates through magnetic dipole torques and GW emission, is given by: 

dE 9W _ Kiy3{1 + with y e [Q _ 2/ p o] (13) 



dv n 2 I 



zz 



where 



5c 2 R e B 2 sin 2 a 1 } 

In this expression R is the radius of the star, e = (I xx — I yy )/I zz the ellipticity, Iij the 
principal moment of inertia , B the magnetic field and a the angle between the rotation 
and the dipole axis. 

In the original scenario of [34, 35], super-strong crustal magnetic fields {B ~ 
10 14 — 10 16 G) can be formed by dynamo action in proto neutron stars with very small 
rotational periods, larger than the break up limit around 0.5 — 1 ms, but smaller than the 
convective overturn at 3 ms. For these highly magnetized neutron stars, the distortion 
induced by the magnetic torque, becomes significant, overwhelming the deformation due 
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to the fast rotation. In the case when the internal magnetic field is purely poloidal and 
matches to the dipolar field B in the exterior, is given by [37]: 

p R 8 B 2 sin 2 a 

es = P 4GIl (15) 

According to the numerical simulations of [37], the distortion parameter /3, which 
depends on both the equation of state (EOS) and the magnetic field geometry, 
can range between 1 — 10 for a non-superconducting interior to 100 — 1000 for a 
type I superconductor and even take values larger than 1000 — 10000 for a type 
II superconductor with counter rotating electric currents. Assuming R — 10 km, 
I zz = 1.4 x 10 45 g cm 2 , P = 1 ms, (3 = 100 and the average value given by the 
observation of SGR and AXP B = 10 15 G , we obtain £ B ~2x 10" 4 and K ~ 3 x 10 37 
erg Hz -3 . In this case, the GW emission becomes negligible compared to the magnetic 
torque and eq. 13 simplifies to 

^ ~ K v 3 (16) 

Considering that magnetars represent about 10% of newborn NS, in agreement with the 
estimates of [36] and population synthesis of [17], we obtain for the mass fraction of the 
progenitors A = 0.1 where the mass fraction of neutron star progenitors is given by: 

/■40M© 

\ NS = / £(m)dm = 9 x lO^M" 1 (17) 

J8M 

where £(m) is the modified A Salpeter IMF from eq. 7 and where we have assumed that 
NS progenitors have masses larger than 8 M [28] and that stars with masses larger 
than 40 M give rise to black holes. Evolution of such massive stars being very fast, we 
can replace the previous expressions in eq. 10, and obtain: 



^^xlO-^f^'^^l 1 ^ (18) 

The energy density increases as z/ 4 at low frequencies and reaches a maximum of 
fl p ~lx 10~ 10 around 1100 Hz. 

When the product f3B is large enough (K » 7r 2 iV~ 2 ), GW emission becomes 
the most important process. In the saturation regime where the spindown is purely 
gravitational, eq. 13 simplifies to: 

^ ~ n 2 Iu (19) 
dv 

and the energy density increases as v 2 at low frequencies and reaches a maximum of 
VL gw ~ 1.3 x 10~ 8 around 1600 Hz (Fig. 1). It has been suggested that the internal 
magnetic field could be dominated by the toroidal component [39, 40]. In this case 

e B = 1.6 x 10~ 4 < B 2 t m > (20) 

where < B 2 16 > is the mean value of the internal toroidal component in unit of 10 16 G, 
and whether the rotational energy is dissipated due to dipole or GW emission depends 
on the ratio B 2 /B. 
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Figure 1. energy density of the background produced by tri-axial rotating magnetars, 
as a function of the observed frequency, for the canonical model with a poloidal internal 
magnetic field and ms initial rotational period. The black continuous line corresponds 
to a model with B = 10 15 G and g = 100 (magnetic spindown regime), the grey 
continuous line to the purely gravitational spindown regime and the black dashed line 
to an intermediate case with B = 10 16 G and g = 1000. 



4-2. Double neutron star coalescence 

The merger of two neutron stars, two black holes or a neutron star and a black hole 
are among the most important sources of gravitational waves, due to the huge energy 
released in the process. In particular, double neutron stars (DNSs) may radiate about 
10 53 erg in the last seconds of their inspiral trajectory. In a recent work, [15, 16] 
used Monte Carlo simulations to calculate the contribution of DNSs to the stochastic 
background, in the frequency band of ground based interferometers, which corresponds 
to the last ~ 1000 s before the last stable orbit, when more than 96% of the gravitational 
energy is released. At that time, the system has been circularized through GW emission 
and the spectral energy density is given in the quadrupolar approximation by: 

gw = K b v- l ' z with v e [10 - v LSO ] (21) 



dv 



where 



(Gtt) 2 / 3 mi m 2 

For double neutron stars with masses m\ = m 2 = 1.4 M , one obtains K b = 5.2 x 10 50 
erg Hz -2 / 3 and the gravitational frequency at the last stable orbit is assumed to be 
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v LSO = 1.57 kHz [41]. 

The merging occurring long after the formation of the system of massive stars, the 
coalescence rate per unit of comoving volume, which will replace R*(z) in eq. 5 and 
eq. 10 is given by: 

R c {z) = f ^^R*(t c (z) -t d )P(t d )dt d (23) 

J L + Zf 

where Zf is the redshift of formation of the progenitors and t c (z) the cosmic time at 
the redshift of coalescence z. The probability distribution of the delay time is usually 
parameterized by [42, 43, 44]: 

P d (t d ) oc i with t d > r (24) 

where the minimal delay time r ~ 20Myr corresponds roughly to the time it takes 
for massive binaries to evolve into two neutron stars [45]. The mass fraction converted 
into the progenitors is given by the product A& = ftNsfb^NS, where /3ns is the fraction 
of binaries which remains bounded after the second supernova event, /& the fraction of 
massive binaries formed among all stars and X^s the mass fraction of NS progenitors. 
In our reference model we assume (3ns h = 0.03 [42]. The merging rate in our galaxy at 
time t is given by: 

rmw(t) = At / p mw (t - t d )P d (t d )dt d (25) 

J To 

where p mw (t) = n^R mw (t) is the galactic star formation rate in M yr _1 . Using the 
model of [46] for the star formation history in the disk, namely 

R mw {t) = 0.056^e"A (M yr- 1 Mpc- 3 ) (26) 

where t is in Gyr, and assuming a density of milky-way equivalent galaxies of n mw = 0.01 
Mpc~ 3 , we find an actual rate of r° lw = r mw (13.5) ~ 3 x 10~ 5 yr _1 , in agreement with 
the expectations derived from statistical studies and population synthesis [47, 42, 48]. 
We find that sources at redshifts z > 0.5 contribute to a truly continuous stochastic 
background, while sources at redshifts 0.25 < z < 0.5 are responsible for a popcorn 
noise, with duty cycles of 1 and 0.1 respectively. The energy density 

A _ o a w m-10,.2/3 f*™* R c( z ) 

1 L n 



ro) = 8, Xl O-^3 r _^i_ d2 (27) 



reaches a maximum of Q gw ~ 7 x 10 -10 around 500 Hz for the continuous contribution 
and of Vt gw ~ 9 x 10 -10 around 550 Hz for the popcorn background (Fig. 2). The 
total background, including the nearest sources down to z ~ is slightly higher, with a 
maximum of Vt gw ~ 1.2 x 10 -9 at 600 Hz. 



5. Constraints with the next generations of detectors 



We now discuss the constraints one can expect to put on some simple models, with the 
next generations of terrestrial detectors. 
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Figure 2. energy density of the continuous background produced by DNS coalescences 
at z > 0.5 (continuous line) and of the popcorn contribution corresponding to sources 
between z = 0.25 — 0.5 (dashed line). The signal from the whole population down to 
z = is also plotted for comparison (dot-dashed line). 



( iv } -5o^ T 7o df fwmm (28) 



5.1. detection 

The optimal strategy to search for a stochastic background, which can be confounded 
with the intrinsic noise background of the instrument, is to cross correlate measurements 
of multiple detectors. The optimized S/N ratio for an integration time T is given by 
[29]: 

where P\{f) and P2C/O are the power spectral noise densities of the two detectors and 
7 is the normalized overlap reduction function, characterizing the loss of sensitivity due 
to the separation and the relative orientation of the detectors. We use the advanced 
LIGO sensitivity as an example of the second generation, the Einstein Telescope as an 
example of the third generation, and the LIGO Hanford/Livingston pair (Hl-Ll) as an 
example of separated detectors. For separated detectors the signal to noise ratio {S/N) 
is calculated between 10 — 150 Hz and for co-aligned and coincident detectors, between 
10 — 500 Hz. We assume an observation time of T = 3 yr and detection threshold of 
S/N = 1 and S/N = 5. 



Astrophysical Sources of Stochastic Gravitational- Wave Background 



9 



5.2. magnetars 

We investigate the constraints on the mean value of the effective external magnetic 
field B e ff and the deformation parameter (3, for the canonical model with ms initial 
periods (dynamo process). Fig. 3 shows the detection limit in the plane B e ff — f3 for 
advanced and third generation co-aligned and coincident detectors and for thresholds 
S/N — 1 and 5. At least with the third generation, we should be able to exclude the 
highest possible values of the magnetic field (B c g > 10 16 G) and models with type II 
superconducting interior or with counter rotating electric currents. We could also put 
constraints on the internal toroidal component of the magnetic field B t . If there is no 
detection, we should be able to confirm or rule out models for which the spindown is 
dominated by GW emission, for instance with B t » B [40], which give S/N ~ 2 with 
advanced detectors and S/N ~ 23 with the third generation. With the first generation 
of detectors, unfortunately, it won't be possible to say anything, since we are well below 
the detection threshold, even for the most optimistic case when the spindown of the star 
is due to GW emission only. 

5.3. Double neutron star coalescence 

We use the model of section 4 to investigate the constraints on the fractions ft, and 
(3 ns , or equivalently on the galactic coalescence rate R m w (Table 1), given by statistical 
studies and source modeling in the range — 10~ 6 yr _1 , more likely around 10 -5 
yr _1 [47, 42, 48]. With advanced coaligned and coincident detectors, we should be able 
to exclude the most optimistic predictions. It would become extremely interesting with 
the third generation, since we expect to reach the actual theoretical expectations with 
separated detectors and get very close to the lower limit with co-aligned and coincident 
detectors. If there is no detection with the third generation, we may have to review our 
models of binary evolution. 

6. Conclusions 

In this article, we reviewed the spectral and statistical properties of astrophysical 
backgrounds and presented two promising models for ground based detectors, as well 
as the constraints one can expect to put on the source parameters with the next 
generation of detectors. We showed that the energy density of astrophysical backgrounds 
is characterized by an increase at low frequencies (usually a power law) , a maximum and 
a cutoff, at the maximum emission frequency in the source frame. For most of our models 
the peak occurs at kHz frequencies, where the sensitivity of pairs of separated detectors 
drops significantly. The best strategy may be to use co-aligned and coincident detectors 
such as the two LIGO Hanford interferometers; the issue is the presence of correlated 
noise, but new techniques are under development in the LIGO collaboration which 
should allow to use this pair in the near future [51]. On the other hand, the stochastic 
background may result in a popcorn noise or may be anisotropic at close redshifts, but 
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Figure 3. magnetars: constraints on the dipole external magnetic field B c j$ and 
the distortion parameter (3 with co-aligned and coincident detectors and integration 
time T = 3 yr. The black continuous line corresponds to advanced detectors and 
a detection threshold S/N = 1, the light grey lines to 3rd generation detectors and 
S/N — 1 (continuous) and S/N — 5 (dashed). The grey continuous line shows the limit 
above which the spindown is almost purely gravitational and corresponds to S/N ~ 2 
with advanced detectors and S/N <~ 23 with the third generation. The average values 
of B e g for observed AXPs and SGRs as well as the magnetar limit, are also plotted 
for comparison. 





(S/Ny l f b p ns 


(S/N)- l R mw 


yr 1 


R mw yr- 1 (S/N = 5) 


inital (HL) 


138.3 


1.3 




6.6 


initial (c-co) 


1.55 


0.015 




0.074 


Ad (HL) 


0.047 


4.5 x 10" 


4 


0.0022 


Ad (c-co) 


0.0025 


2.4 x 10- 


5 


1.2 x 10" 4 


3rd (HL) 


4.5 x 10" 4 


4.3 x 10- 


6 


2.2 x 10- 5 


3rd (c-co) 


1.6 x 10" 4 


1.5 x 10- 


6 


7.5 x 10- 6 



Table 1. Constraints on the product fb(3 ns and on the galactic coalescence rate R mw of 
double neutron stars for different generation of detectors, for an integration time T = 3 
yr and detection thresholds S/N = 1 (column 1 and 2) and S/N = 5 (third column). 
HL indicates a pair of detectors separated such as the LIGO Hanford-Livingston pair, 
and c-co a pair of co-aligned and coincident detectors. One can obtain the constraints 
for any detection threshold from columns 1 and 2 by multiplying the values given in 
column 1 and 2 (S/N = 1) by the S/N. Column 3 corresponds to S/N = 5) 
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again adequate detection strategies are also under investigation [52, 53]. We showed 
that we could already obtain interesting astrophysical results on simple models with 
advanced and third generation detectors, such as upper limits on the magnetic field or 
the ellipticity of magnetars or the coalescence rate of compact binaries. On the other 
hand we expect to be able to rule out or confirm some extreme models such as pure GW 
spindown in magnetars. Work is currently in progress to investigate a larger range of 
models and parameters and it will be reported in a future paper. What is particularly 
interesting with stochastic backgrounds, is that we can put constraints on the mean 
value of the population, and not just on some particular sources that may be in the tail 

of the probability distributions of the ellipticity or the magnetic field. 
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